Measurement of angle (3 with time-dependent CP asymmetry in B° — > K + K K° decays 
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We present recent results on CP-violation, and the determination of CKM angle /3, with the decay 
B° -» K + K~K°, with BABAR and Belle detectors. 
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I. INTRODUCTION 

In the Standard Model (SM) of particle physics, 
the phase of the Cabibbo-Kobayashi-Maskawa (CKM) 
quark-mixing matrix [l], 0| is the only source of CP vi- 
olation in the quark sector. Due to the interference be- 
tween mixing and decay, this phase can be observed in 
measurements of time-dependent CP asymmetries of B° 
mesons. In the SM, CP asymmetries in b — > sss decays, 
such as B° — ► K + K~K °, are expected to be nearly equal 
to those observed in tree-dominated b — ► qs decays §■ 
However, because in the SM the former are dominated 
by loop amplitudes, new particles in those loops poten- 
tially introduce new physics at the same order as the SM 
process. Within the SM, deviations from the expected 
CP asymmetries in B Q — ► K + K~K° decays depend on 
the Dalitz plot position, but are expected to be small 
and positive 0]. In particular, for the decay B° — > (j)K° 
they are expected to be less than 4%. BABAR extracts 
the time-dependent CP-violation parameters by taking 
into account different amplitudes and phases across the 
B° and B° Dalitz plots, while Belle measures it sepa- 
rately for B° -> <pK and the rest of K+K~K° events, 
neglecting interference between intermediate states. 

The analyses presented here are based on 347 (535) 
million BB pairs collected with the BABAR (Belle) de- 
tector at the SLAC PEP-II (KEKB) e+e~ asymmetric- 
energy collider. Data are collected on the T(4S) res- 
onance, while a fraction of about 10% is collected at 
approximative^ 40 MeV below the T(AS) resonance, 
and it is used to study the background arising from 
e + e~ — > qq(q = u,d,s,c) continuum events. The BABAR 
and Belle detectors are described in detail elsewhere 



II. EVENT RECONSTRUCTION 

Events are fully reconstructed combining tracks and 
neutral clusters in the detector to form B° — > K + K~K °, 
with a K° reconstructed as K° — > tt + tt~ (B° {+ _ ) ) (BABAR 
and Belle) and K° s -> 7r°7r° (B° (00) ), or K° L (B° (L) ) 
(BABAR only) . In order to select B candidates we use a set 
of two kinematic variables: the beam-energy-substituted 



mass m E s = \/( s / 2 + P» -Pb) 2 /E? 
and the energy difference AE = 
(Ei,pi) is the four-vector of the initial e + e~ system, y/s 



f p 2 B (M bc for Belle), 
E% - y/s/2. Here, 
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is the center-of-mass energy, is the reconstructed mo- 
mentum of the B° candidate, and E* B is its energy calcu- 
lated in the e + e~ rest frame. For signal decays, the raES 
distribution peaks near the B° mass with a resolution of 
about 2.5 MeV/c 2 , and the AE distribution peaks near 
zero with a resolution of 10 — 50 MeV, depending on the 
final state. For decays with a K®, K® momentum is not 
measured, but evaluated by constraining the B° mass to 
the nominal value [6( In this case only AE is used, and it 
has a resolution of about 3 MeV. The main background 
comes from random combinations of particles produced 
in continuum e + e~ — > qq. In the CM frame, these events 
have a jet-like structure, while B decays have a nearly 
isotropic topology. We parameterize this difference us- 
ing several variables, providing additional discrimination 
between signal and background. Another source of back- 
ground comes from decays of B mesons which mimic the 
signal. This background is typically more difficult to sup- 
press. The contribution of these decays is estimated from 
Monte Carlo simulations. 

For each fully reconstructed B° meson (Bcp), we use 
the remaining particles in the event to reconstruct the de- 
cay vertex of the other B meson (-Btag) and identify its 
flavor <?tag- A multivariate tagging algorithm determines 
the flavor of the B tag meson in BABAR data and classi- 
fies it in one of seven mutually exclusive tagging cate- 
gories depending on the presence of prompt leptons, one 
or more charged kaons and pions Q ■ The performances of 
this algorithm are measured with a data sample of fully 
reconstructed B° decays into flavor eigenstates (-Bg av ): 
B° — > £>(*)~7r + / p + /a^ . The effective tagging efficiency 
is Q = X; c e c (l~2w c ) 2 = 0.304±0.003 for BABAR (similar 
for Belle), where e c (w c ) is the efficiency (mistag proba- 
bility) for events tagged in category c. For Belle, the tag- 
ging algorithm returns qt ag and the tag quality r, which 
varies from r = for no flavor discrimination to r = 1 
for unambiguous flavor assignment. Events with r < 0.1 
are discarded for the CP-asymmetry measurement, and 
the others are sorted into six intervals. The difference 
At = tcp — h ag of the proper decay times of the Bcp 
and P>tag mesons is calculated from the measured dis- 
tance between the reconstructed decay vertices and the 
boost (/3 7 = 0.56(0.465) for BABAR (Belle)) of the T{AS). 
The error on At, a At, is also estimated for each event. 
Events are accepted if the calculated At uncertainty is 
less than 2.5 ps and |At| < 20 ps. The fraction of events 
which satisfy these requirements is 95%. 

The BABAR analysis strategy is to perform a maxi- 
mum likelihood fit to the selected K + K~K° events with 
a likelihood function C, which uses as probability den- 
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sity function (PDF) for each event, C = ^(m^^AE) • 
Plow • V DP {m K + K - , cos6 H , At, q tag ) x Tl(At, <j At ) where 
rii is the yield for each category (i = signal, continuum 
background, and BB backgrounds), and 1Z is a At reso- 
lution function with parameters determined in the -Bfl av 
data sample. Plow is a supplementary PDF used only 
in the fits to the region with m K + K - < 1.1 GeV/c 2 dis- 
cussed below. It depends on the event shape variables 
and, for B° (I/) only, the missing momentum of the event. 
This PDF accounts for the fact that signal decays have a 
missing momentum consistent with the reconstructed 
direction, while for background events it is more isotrop- 
ically distributed. We characterize events on the Dalitz 
plot in terms of the invariant K + K~ mass, m K + K ~ , and 
the cosine of the helicity angle between the K + and the 
K° in the CM frame of the K + K~ system, cobOh- 

In the Belle approach the likelihood for the event selec- 
tion is the same, without Vdp- A loose requirement on 
the likelihood ratio H s /b = £sig/{£sig + £-bkg) is applied, 
and a maximum-likelihood fit to observed At distribution 
is performed to the selected events. 



III. ANALYSIS OF DALITZ PLOT 

Accounting for the experimental efficiency e, the Dalitz 
plot PDF for signal events is 

Vdp = dT ■ e(m K + K -,cos9 H ) ■ \J(m K+K -)\, (1) 

where the final term is the Jacobian of the transformation 
for our choice of Dalitz plot coordinates. 

The time- and flavor-dependent decay rate over the 
Dalitz plot is 



dT 



,-\At\/T 

At 



\A\* + \A\* (2) 
Vcp <?ta g 2Im(AA*e~ 2if3 ) sinAm d Ai 
q (\A\ 2 - \A\ 2 ) cos ArridAt 



where tjcp = 1 (-1) for B° (+ -)> B V> (5%)- T and 



Arrid are the lifetime and mixing frequency of the B° 
meson, respectively [|| . The CKM angle {3 enters through 
B°-B° mixing. Actual world average is /3 ~ 0.38 [f|. We 
define the amplitude A (A) for B° (B°) decay as a sum 
of isobar amplitudes, 



A(A) = ^c r (l±6 r )e 4 



r±S r ) 



/(/), 



(3) 



where the parameters c r and <p r are the magnitude and 
phase of the amplitude of component r, and we allow 
for different isobar coefficients for B° and B° decays 
through the asymmetry parameters b r and 5 r . Our 
model includes the vector meson 0(1020). We include 
also decays into intermediate scalar mesons: /o(980), 
A o (1550), and Xco- The angular distribution is con- 
stant for scalar decays, whereas for vector decays is 



Z r = — 4<f ■ p, where q is the momentum of the res- 
onant daughter, and p is the momentum of the third 
particle in the resonance frame. We describe the line- 
shape for the 0(1020), A (1550), and Xco using the rel- 
ativistic Breit-Wigner function [8]. For the 0(1020) and 
Xco parameters, we use average measurements @. For 
the Ao(1550) resonance, we use parameters from the our 
analysis of the B+ -> K+R-R+ decay @. The / (980) 
resonance is described with the coupled-channel (Flatte) 
function Q, with the coupling strengths for the KK 
and 7T7T channels taken as g n = 0.165 ± 0.018 GeV/c 2 , 
9k/9tt = 4.21 ± 0.33, and the resonance pole mass 
m r = 0.965 ± 0.010 GeV/c 2 0. In addition to resonant 
decays, we include three non-resonant amplitudes. The 
existing theoretical models, which consider contributions 
from contact term or higher-resonance tails [l(| [HI, [12] do 
not reproduce well the features observed in data. There- 
fore we adopt a phenomenological parameterization [l3T ] 
and describe the non-resonant terms as an exponential 
decay: 

A NR = c ^ e ~ Qm ' J • (1 + b NB.) ■ e^+^> (4) 

and similarly for Anr- 

We compute the CP-asymmetry parameters for com- 
ponent r from the asymmetries in amplitudes (b r ) and 
phases (5 r ) given in Eq. ([3]). The rate asymmetry is 



.4, 



CPr 



| Ar j | Ar 



-2b r 



i + br 



(5) 



|A| 2 + |A| 2 

and (3 e ff — (3 + 5 r is the phase asymmetry. 
The fraction for resonance r is computed 

Jdcos6 H dm K+K - ■ \J\ ■ {\A r \ 2 + |A| 2 ) f . 
Jdcos8 H dm K+K - ■ \J\ -(\A\ 2 + \A\ 2 ) ' U 

The sum of the fractions can be larger than one due to 
negative interference in the scalar sector. 

The fit to BaBAR data returns 879±36 B° (+ _ ) , 138±17 
B° (00) and 499 ± 52 B° (L) signal candidates. The isobar 
amplitudes, phases and fractions are listed in Table [U 
Signal weighted distribution for the Dalitz plot projec- 
tions in the entire phase space and in a reduced region 
m K+K- < 1-1 GeV/c 2 , where we extract separate CP 
asymmetry parameters, are shown in Fig. [TJ 

The fit to Belle data returns 840 ± 34 signal B° (+ _ } 
candidates, after vetoing the with the requirement 
\m K +K- ~ m 4>\ > 15 MeV/c 2 . 

As a cross-check of the Dalitz model extracted by the 
fit, we compute angular moments and extract strengths 
of the partial waves in m K + K - bins using the B° (+ _ ) 
sample. In this approach we only assume that the two 
lowest partial waves are present. We verified this as- 
sumption determining that the higher angular moments 
((P 3 _ 5 )) are consistent with zero. In our model, the 
P-wave contribution comes from 0(1020)^° decays and 
from non-resonant events with K + K® and K~K® mass 
dependence. We find that the total fraction of P-wave in 
the entire Dalitz plot is f p = 0.29 ± 0.03(stat). 
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TABLE I: Isobar amplitudes, phases, and fractions from the 
fit to BABAR, data. The fraction for non-resonant amplitude 
is given for the combination of the three contributions. 



Decay 


Amplitude c r 


Phase 4> r 


Fraction Tr (%) 


0(lO2O)if u 


0.0098 ± 0.0016 


-0.11 


±0.31 


12.9 ± 1.3 


fo(980)K° 


0.528 ±0.063 


-0.33 


±0.26 


22.3 ± 8.9 


X (1550)A"° 


0.130 ± 0.025 


-0.54 


±0.24 


4.1 ± 1.8 


NR (K+K-) 


1 (fixed) 





(fixed) 




(K+K°) 


0.38 ±0.11 


2.01 


±0.28 


91 ± 19 


(K~K°) 


0.38 ±0.16 


-1.19 


±0.37 




XcoK u 


0.0343 ± 0.0067 


1.29 


±0.41 


2.84 ± 0.77 


D+K~ 


1.18 ±0.24 






3.18 ±0.89 


D+K- 


0.85 ±0.20 






1.72 ± 0.65 



> 250 r 






FIG. 1: Distributions of the Dalitz plot variables (left) 
m K + K - and (right) cos 8h for signal events (points) com- 
pared with the fit PDF for candidates. Top distri- 
butions are for the entire phase space, bottom for a reduced 
region m K+K - < 1.1 GeV/c 2 . 



IV. CP ASYMMETRY 

In the Belle "quasi-two-body" approach, the time de- 
pendent decay rate in Eq. [3] is simplified because the in- 
terference effects are neglected: 



/±(At) 



; -|At|/r 



(7) 



[ 1 ± S sin (Am d At) =f C cos (Am d At) } 



The parameters C and S describe the amount of CP vi- 
olation in decay and in the interference between decay 
with and without mixing, respectively. The SM expec- 
tations for B° -> K+K-K° s are S = -(2/+ - l)sin2/3, 
where f+ is the CP-even fraction. Using isospin rela- 
tions, /+ has been measured on 357 fb _1 data sample, 
and gives /+ = 0.93 ± 0.09(stat) ± 0.05(syst) [H. This 
result is confirmed by partial wave analysis performed on 
BABAR data [3]. 

In the BABAR analysis, the time-dependent fit to the 
Dalitz plot allows to extract (3 e ff removing the trigono- 
metrical ambiguity (3 e ff — > 7r/2 — f3 e ff- In this analysis 



TABLE II: Time-dependent CP-asymmetries Acp and j3 e ff for 
B° -> K + K~ K° in m K+K - < 1.1 GeV/c 2 and in the whole 
phase space (BABAR), and C and S for B° -» K + K~ K% in 
m K + K - > 1.034 GeV/c 2 (Belle). The first error is statistic, 
the second systematic. The third error in S is systematic 
effect due to the CP-content uncertainty. 



Decay 



A C p(<t>K l >) 
A C p(foK°) 

MM ) 



CP asymmetry 



-0.18 ±0.20 ±0.10 
0.06 ±0.16 ±0.05 
0.45 ±0.28 ±0.10 
0.18 ±0.19 ±0.04 



Acp(K + K~ K u ) 
f3 eff (K+K-K°) 
C(K + K'Kl) 
sm(2/3 e g)(K + K-K°) 



-0.034 ± 0.079 ± 0.025 
0.361 ± 0.079 ± 0.037 



0.09 ±0.10 ±0.05 
0.68 ± 0.15 ±0.03!£; 2 3 



the reflection is suppressed from the interference between 
CP-even and CP-odd decays that give rise to a cos(2/3 eJ y) 
term in EqjJH in addition to the sin(2/3 eJ y) terms that 
come from the interference decays with and without mix- 
ing. In this case we measure an average (3 and Aqp for 
the full K + K~K° phase space. 

We measure also CP-asymmetry parameters in the re- 
gion m K + K - < 1.1 GeV/c 2 , where the model depen- 
dence from the rest of the Dalitz plot is highly reduced. 
In this region, we fit the CP time-dependent asymme- 
tries for the </>(1020) and /o(980) components, while we 
fix the ones for the low-K + K ~ tail of the non-resonant 
decays to the SM expectation. The Dalitz plot model 
is fixed to the one measured in the full phase space and 
reported in Table HI with the exception of the <^>(1020) 
isobar coefficients, which are fitted simultaneously with 
CP-asymmetry parameters. In this reduced region of the 
phase space we find 252 ± 19, 35 ± 9, 195 ± 33 signal 
events for B° (+ _ ) , B° (00) and B° (L} respectively. The re- 
sults on the CP asymmetries are shown in Table [TTJ The 
sources of systematic uncertainties are briefly described 
below. 

The significance of the nominal result for f3 e ff in the 
entire Dalitz plot, compared to the trigonometrical re- 
flection is of 4.6cr. The significance of the CP-violation is 
4.5c. In Fig. [2] the distributions of At for £?°-tagged 
and -B°-tagged events, and the asymmetry A(At) = 



(N B o — Ng )/(N B o + Ngo), for background subtracted 
BABAR data are shown. In Fig. [3] the asymmetries for 
good-tagged (r > 0.5) K + K~K® events in Belle data 
are shown. 

Systematic effects are associated to parameterization 
of the signal PDFs, possible fit bias, fixed At resolu- 
tion parameters, B° lifetime, B°-B° mixing and flavor 
tagging parameters. Smaller errors due to beam-spot 
position uncertainty, detector alignment, and the boost 
correction are based on studies done in charmonium de- 
cays. In BABAR analysis, an uncertainty is assigned to 
the resonant and non-resonant line-shapes. We try sev- 
eral alternative non-resonant models which omit some 
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FIG. 2: (top) At distributions and (bottom) asymmetries for 
events for (left) 1.0045 < m K+K - < 1.0345 GeV/c 2 
and (right) the whole Dalitz plot for BABAR data. For the 
At distributions, B°- (B -) tagged signal-weighted events are 
shown as filled (open) circles, with the PDF projection in solid 
blue (dashed red). 
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FIG. 3: Asymmetry for good-tagged events (r > 0.5) for 
B° -> K + K~K° for Belle data. The solid curve shows the 
PDF projection for the result of the unbinned maximum- 
likelihood fit. The dashed curve shows the SM expectation 
with the measurement of CP-violation parameters for the 
B° -> J/ipK° decays. 

of the dependencies on K + K° and K~ K° masses (see 
Eq. We also study the effect of the uncertainty of 
the shape parameter a on the CP parameters. The non- 
resonant events contribute to the background under the 



(j> but their shape is determined from the high-mass re- 
gion. We therefore omit the non-resonant terms and take 
the difference from the reference fit as a systematic er- 
ror for 4>K° and f o (§&0)K° CP-asymmetries. In Belle 
measurement, by far the largest systematic uncertainty 
on CP parameters is associated to the knowledge of the 
CP-even fraction /+. 

V. CONCLUSIONS 

We have measured the time-dependent CP- 
asymmetries in B° — > K + K~K° decays, with a 
simultaneous analysis of the Dalitz plot distribu- 
tion of the intermediate states (BABAR) or with a 
"quasi-two-body" approach (Belle). The measured 
value of CP- asymmetries in the entire Dalitz plot is 
(3 eff = 0.361 ± 0.079 ± 0.037, which is consistent with 
the SM expectations (f3 ~ 0.38). The trigonometrical 
ambiguity in (3 e ff is removed at 4.6<r. This is the first 
such measurement in penguin modes. Additionally, we 
extracted the CP-asymmetry parameters for B° — > <f>K° 
and B° -> f (890)K°, to be f3 eff = 0.06 ± 0.16 ± 0.05 
and l3 eB = 0.18 ± 0.19 ± 0.04, respectively. Therefore 
we do not observe significant deviations from the SM 
predictions. 

In Belle measurement, the measured sin(2/3 eJ y) = 
0.68 ± 0.15 ± O.Q3±g;?| for B° -> K+K~K° decays with 
exclusion of <f>K® events is also consistent with SM pre- 
dictions. 
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